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The annihilation of neutraUno dark matter may result in observable signals in different 
wavelength. In the present paper we will discuss the effect of neutralino annihilation 
in the halo of our Galaxy and in its center. According to high resolution cold dark 
matter simulations, large virialized halos are formed through the constant merging of 
smaller halos appeared at previous times. At each epoch, dark matter halos have then a 
clumpy component which is made of these merging subhalos. The annihilation of dark 
matter in these clumps, always present in the halo of our Galaxy, may be responsible 
for appreciable fluxes of 7-rays, potentially detectable. We find that, depending on the 
fundamental parameters of the clump density profile and on the distribution of clumps 
in the Galactic halo, the contribution to the diffuse 7-ray background from clumps could 
be used to obtain constraints on the neutralino properties such as mass and annihilation 
cross section. On the other hand the annihilation of neutralino dark matter in the galactic 
center may result in radio signals. At the galactic center, infact, the accretion flow 
onto the central black hole sustains strong magnetic flelds that can induce synchrotron 
emission, in the radio wavelength, by electrons and positrons generated in neutralino 
annihilations during advection onto the black hole. We flnd that the observed emission 
from the galactic center is consistent with neutralinos following a Navarro Frenk and 
White density proflle at the galactic center while it is inconsistent with the presence of 
a spike density profile, supposed to be generated by the formation history of the central 
black hole. 



1. Introduction 

Most of the matter in tlie universe has yet to be observed in any frequency band, 
thus tlie name, dark matter (DM). The evidence for the predominance of dark 
over visible matter comes mainly from the gravitational effects of the dark matter 
component. However, gravitational studies have been unable to shed light on the 
nature of the dark matter. Big bang nucleosynthesis constrains most of the dark 
matter to be of non-baryonic origin. This has encouraged the study of plausible 
new particle candidates for the dark matter. 

Weakly interacting massive particles (WIMPs) are natural candidates for the 
dark matter. Particles with masses around ^ 100 GeV that interact only weakly 
have freeze-out densities in the required range of densities. In addition, particle 
physics models that invoke supersymmetry generate a number of plausible WIMPs. 
In the supersymmetric extensions of the standard model, the lightest supersymmet- 
ric particle may be stable due to conservation of R-parity enabling their survival to 
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the present. In addition to massive and weakly interacting, dark matter particles 
are expected to be neutral. A class of neutral lightest supersymmetric particles is 
represented by a combination of gauginos and higgsinos, named the neutralino often 
represented by x- 

Given the requirements of neutralino production in the early universe, it ispos- 
sible to study the phenomenology of such dark matter candidates in detail!^. In 
particular, the annihilation of neutralinos has often been considered a potential 
source of detectable secondaries: high energy particles and electromagnetic radia- 
tion. In this sense, dark matter can be visible through the radiation caused by the 
annihilation secondaries E 3 4 o 6 7 8 9 10 11 12 13 14 15 16 gj^ce the neutralino is 
a Majorana particle, it will self-annihilate at a rate proportional to the square of 
neutralino density. Thus, as was first realized in , the highest density dark matter 
regions are the best candidates for indirect searches. In the present paper we will 
review the annihilation signal that could come from the clumped halo component of 
our Galaxy and from its Galactic Center (GC), in these regions infact it is expected 
a large enhancement in the neutralino density with a consequent amplification of 
the annihilation signal. 

Recent advances in Cold Dark Matter (CDM) simulations have shown that the 
large scale structure of the Universe can be explained in terms of a hierarchical 
scenario in which large halos of dark matter are generated by the continuous merging 
of smaller halos this picture a dark halo is the superposition of a smooth 

component, characterized by a typical scale comparable with the virial radius of 
the forming structure, and a clumped structure made of thousand of small scale 
halos. CDM simulations also show that most halos are well described by a density 
distribution with cusps at the center of each halo. The exact shape of the central 
cusp is still a matter of debate. Most recent simu lations favor profiles with density 
cusps varying from the Moore et al. profile where pDAiir 0) ^ ^-i-S ^ 
Navarro, Frenk, and White (NFW) profile "221 where pdm{t — > 0) ~ r^^ . 

The debate is exacerbated by observations of galaxy rotation curves that seem to 
provide no evidence of central cusps However, the survival of cusps in galactic 
centers is highly dependent on the galaxy's mer ger history in particular on the 
history of formation of galactic center black holes The central regions of small 
mass dark matter halos (DM clumps) are less affected by the dynamics of baryonic 
matter and less likely to have black hole (BH) mergers at their centers. Thus, a 
cuspy profile may well describe the density of DM clumps. 

Another important piece of information is embedded in the spatial distribution 
and survival history of clumps on their way to the central part of the host galactic 
halo. Much physics is involved in the description of the structure of a DM clump 
moving in a larger DM halo, and different recipes are possible. These unknowns 
have forced us to consider two different scenarios, that for simplicity we call type I 
and type II scenario. 

In the type I scenario, the clump position in the host Galaxy determines its 
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external radius. In the case of NFW and Moore profiles, the location of the core 
is then determined by assuming a fixed fraction of this external radius. In the 
type II scenario, a recipe is taken from the literature for the so-called concentration 
parameter, defined as the ratio of the core and virial radii of a clump. The recipe 
allows one to determine the properties of a clump for a given mass. We show that 
in this scenario the clumps are gradually destroyed on their way to the center of 
the Galaxy, so that an inner part, depleted of its clumpy structure is formed. This 

no 

seems compatible with numerical simulations . The effects of these two scenarios 
on the 7-ray emission from the annihilation of CDM particles are dramatic: in the 
type I scenario, the highly concentrated clumps that are implied produce strong 
7-ray emission while in the type II scenario the low concentration clumps imply 
7-ray fluxes several orders of magnitude smaller than in the previous case. 

Another important piece of information could come from the annihilation signal 
at the GC. Infact, apart from the clumped halo component, the GC regk)n may 
potentially be so dense that all neutralino models would be ruled out This 
strong constrain arises in models where the super-massive BH at the galactic center 
(GC) induces a strong dark matter density peak called the spike. The existence of 
suc h a spike is strongly dependent on the formation history of the galactic center 
BJjES If the BH is formed adiabatically, a spike would be present while a history 
of major mergers would not allow the survival of a spike. 

In contrast to the uncertain presence of a central spike in the dark matter dis- 
tribution, the central BH is known to induce an accretion flow of baryonic matter 
around its event horizon. The accretion flow carries magnetic fields, possibly am- 
plified to near equipartition values due to the strong compression. The distribution 
of electrons and positrons (hereafter called electrons) produced by neutralino an- 
nihilation at the GC would also be compressed toward the BH radiating through 
synchrotron and inverse Compton scattering off the photon background. 

By considering the injection of electrons, combined with radiative losses and 
adiabatic compression, we find the equilibrium spatial and spectral electron dis- 
tribution and derive the expected radiation signal. We find that the synchrotron 
emission of electrons from neutralino annihilation range from radio and microwave 
energies up to the optical, in the central more magnetized region of the accretion 
flow. At low frequencies, synchrotron self-absorption slightly reduces the amount of 
radiation transmitted outwards. The resulting signal is stronger than the observed 
emission in the 10 to 10^ GHz range for the case of a spiky dark matter profile while 
for a pure NFW profile the emission is below the observed values. 

The paper is structured as follows in the first three paragraphs we will discuss 
the annihilation signal that could come, in the 7-ray frequency range, from the 
clumped halo of the Galaxy, while in the last three paragraphs we will discuss the 
synchrotron signal that could come, in the radio frequency range, from the GC. We 
will conclude in paragraph 7. 
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2. Dark Matter Clumps in the Halo 

The contribution to the diffuse 7-ray background from xx annihilation in the clumpy 
halo depends on the distribution of clumps in the halo and on the density profile 
of these clumps. Our purpose here is to investigate the wide variety of possibilities 
currently allowed by the results of simulations and suggested by some theoretical 
arguments, concerning the density profiles of dark matter clumps. We consider 
three cases: singular isothermal spheres (SIS), Moore et al. profiles, and NFW 
profiles. 

The Moore et al. and the NFW profiles are both the result of fits to different 
high resolution simulations ll^l. Although there is ongoing debate over which profile 
is most accurate, it is presently believed that a realistic descriptions of the dark 
matter distribution in halos will follow a profile in the range defined by the Moore 
et al. and the NFW fits The dark matter density profiles can be written as 
follows: 
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The SIS and Moore et al. clump density profiles are in the form given in eqs. 



Inside rmin, neutralino annihilations 



Q and (0) down to a minimum radius, r 
are faster than the cusp formation rate, so that p(r < rmin) = p(fmin) remains 
constant. To estimate rmin, following 1^, we set the annihilation timescale equal to 
the free-fall timescale, so that 
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where {crv)j^nn is the xX annihilation cross-section, G is the Newton constant, Mj, 
and are respectively the clump mass and the neutralino mass. 

The fundamental parameters of the clump density profile are the density nor- 
malization Po, the clump radius ro and, in the case of Moore et al. and NFW 
profiles, the clump fiducial radius rf. In order to fix these fundamental parameters 
we have considered two different scenarios (type I and II). 

We have modeled the smooth galactic halo density with a NFW density profile 
[Eq. with rf = 27 kpc and po determined from the condition that the dark 

matter density at the Sun's position is pda/ (^0) — 6.5 x lO^^^g/cm . 
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In the type I scenario, the radius of a clump with fixed mass is determined by its 
position in the Galactic halo. More specifically, the radius of the clump is located at 
the radius where the clump density equals the density of the Galactic (smooth) dark 
matter halo at the clump position (namely po). The physical motivation for such a 
choice is to account for the tidal stripping of the external layers of the clump while 
the clump is moving in the potential of the host halo. For the NFW and Moore 
profiles, the fiducial radius rf has been taken as a fixed fraction of tq: rf = O-lrg. 

In the type II scenario the external radius of the clumps is taken to be their 

1/3 

, where p200 is 



3Mc 



virial radius, defined in the usual way: = r„ir — 47rp2oo 
200 times the critical density of the Universe pc = 1.88 x 10~^^/i^ g/cm? (we have 
assumed h = 0.7 everywhere). In this scenario, following^! we have introduced the 
concentration parameter defined as ^ = where r_2 is the radius at which the 
effective logarithmic slope of the profile is —2, set by the equation -j^ p{f) — 0. 
The mass dep endence of the concentration parameter used in our calculations is 
taken fromES. 

The definite trend is that smaller clumps have larger concentration parameter, 
reflecting the fact that they are formed at earlier epochs, when the Universe was 
denser. In general, the concentration parameter has also a dependence on the 
redshift at which the parameter is measured. We are not interested here in such 
dependence, since we only consider what happens at the present time (zero redshift). 
The normalization constant in the clump density profile po is fixed by the total 
clump mass Mc- 

In the case of the NFW density profile rf — r„2, while in the case of the Moore 
et al. density profile rf = r_2/0.63^l In terms of concentration parameters 
^NFW = 0.63^Moore • Usiug the concentration parameter as in for NFW clumps, 
we can estimate the stripping distance, which we define as the typical distance from 
the galactic center where the density of a clump of fixed mass and the density in 
the smooth DM profile are equal at the fiducial radius rf oi the clump. In other 
words, at the stripping distance the layers of the clump outside the fiducial radius 
will have been stripped off. From this estimate it is easy to see that most clumps 
in the inner parts of a host galaxy are stripped off of most of their material, so that 
this region has no clumpy structure. Another way of seeing this phenomenon is 
that the clumps that are able to reach the central part of the Galaxy are effectively 
merged to give the observed smooth dark matter profile. 

The exception to this conclusion may be represented by low mass clumps, which 
are more concentrated (denser cores) and may then penetrate deeper. Numerical 
simulations show the disappearance of large clumps in the centers of galaxy size 
halos, but they cannot resolve the smaller denser clumps that may eventually make 
their way into the core of the galaxy. For simplicity, in our calculations we assume 
that the inner 10 kpc of the Galaxy have no clumps at all. 
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Figure 1. Flux of 7-rays in units of (GeV cm^ s sr) ^ arriving on Earth averaged in all directions 
for = 100 GeV, ((rv)Ann = 3 X lO'^'' cm^/s and Mc^min = 10^ Mq. SIS density profile 
(dotted line), Moore et al. density profile (dashed line) and NFW density profile (continuous 
line). Also shown is the EGRET data on extragalactic diffuse 7-ray background (left panel: first 
scenario; right panel: second scenario). 



3. Gamma-ray emission in Neutralino Annihilation 

In order to determine the 7-ray emission from the DM ch.imps, as aheady pointed 
out in the previous section, the distribution of clumps in the Galaxy is needed. 
The probability distribution function of clumps with a given mass and at a given 
position can be fitted from numerical simulations. In the present paper we follow 
126 adopting a spatial (as a function of the distance d from the galactic center) and 
mass distribution of the clumps reflecting the following expression: 



iV,(d, Me) = N,,0 



Mh 



d 



-3/2 



(5) 



where N^.o is a normalization constant and d^c is the scale radius of the clump 
distribution'^. Simulations find a ~ 1.9 and a halo like that of our Galaxy, with 
Mh ~ 2 X 10^^ Mq, contains about 500 clumps with mass larger than 10^ Mq^I. 
The 7-ray flux per unit solid angle and per unit energy along a fixed line of sight 
in the (0, (f) direction can be computed as 



1 



CMh 



^jiE^,e,(j)) ^ — ds Nc{d{s),M)N^ dM (6) 

where d{s) = ^^ — 2sdQCOs9 + dg is the distance of a generic point on the line 
of sight from the galactic center (with 6 the angle between the direction s and 
the axis Sun-galactic center), M is the clump mass, ^Mh is the maximum allowed 
mass for DM clumps in the Halo (we have used C = 0.01) and is the total 
number of photons emitted per unit time and energy by a DM clump of mass Af . 
This quantity, depending on the scenario chosen for the clump density profile, may 
depend or not on the distance d{s) of the considered point from the galactic center. 



'^In this paper we have assumed dgc = 10 Kpc as in 
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In the first scenario, where the normalization of the clump density is related to 
the smooth Halo density, one has = N^{M, d,Ej), while in the second scenario 

= N^{M, Ej). One should remember that Eq. © is an average over all possible 
realizations of a halo with its clumpy structure. Fluctuations around this value may 
be present due to the accidental proximity of few clumps in the specific realization 
that we happen to experience in our Galaxy. 

We assume that neutralinos mainly annihilate into quark- antiquark pairs, which 
seems co nfirm ed by more detailed calculations carried out in specific supersymmetric 
scenarios El. Actually, it is easy to see that even when neutralinos annihilate into 
pairs of W~^W~ or Z'^Z'^, the end result of the decay chain is dominated by quarks 
and antiquarks, hadronizing mainly into pions (roughly 1/3 neutral pions and 2/3 
charged pions) and their spectral shape is identical to that of direct quark-antiquark 
production. In fact, each W or Z boson would have a Lorentz factor m^/AIz,w 
(if the neutralino mass is large enough to allow the production of a or Z pair). 
In the rest frame of the boson, the maximum energy of the particles generated in 
the decay is Mz^w, so that in the laboratory frame the spectrum has a cutoff at 
m^, as in the case of direct quark production in the neutralino annihilation. The 
spectrum is also left unchanged. 

The number of photons produced with energy Ej in a single XX annihilation 
can be written as follows: 

^'^^-r—,E.PiE.,E/J^ (7) 



dE'y J Ett min ' dE^ 

where P{Et^,E^) — 2{E^ — mj)^^/^ is the probability per unit energy to produce 
a 7-ray with energy E^ out of a pion with energy E^^. For the pion fragmentation 
function we assume the functional form introduced by Hill ^3 

^ = — i|x-3/2(l_x)2 (8) 

with X = ET^/m^, E^^ = and E-^^^i^ = E^ + rn^/AE^. Finally, 

dN,^^ jl-x)^ 

dE^ Am^ J,^ x^^x^ - 772)1/2 ' 

where 77 = m-,r/ni^, and = Ej/m^ + m^rf' / AE.^. 

The neutralino annihilation rate per unit volume is, T^-^{r,M,E^) — 
p\{^) {'^''^) Ann / fn^^, therefore the 7-ray emissivity i-j{r, M, E^) associated to the sin- 
gle clump of mass M is obtained by multiplying Eq. Q by F-^^. The number of 
7-rays produced per unit time and per unit energy in a single DM clump of mass 
M is then 

pro 

N^{M,E^)^ drATTr^j^{r,M,E^) . (10) 
Jo 
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4. Gamma Ray Emission from Clumps 

In this section we present tlie results of our calculations of the 7-ray emission from 
dark matter annihilation in the halo of our Galaxy, including both the smooth 
and clumped components introduced above. We detail the description of these 
results for the two scenarios (type I and II) of spatial distribution of the clumped 
component. We expect the type II scenario to result in a quite weaker signal than 
that obtained in the type I scenario, because the recipe for the type I clumps implies 
much stronger concentration. In the type II scenario the 7-ray signal from clumps 
overcomes the 7-ray flux from the smooth dark matter distribution only in the 
direction of the galactic anticenter and only for SIS and Moore profiles for the dark 
matter distribution inside the clumps. 

On the other hand, in the type II scenario the contribution of the clumped 
component to the diffuse 7-ray flux from the galactic halo is several orders of mag- 
nitude above the smooth halo component, for any of the three density profiles 
considered above. This impressive difference in the predictions is symptomatic of 
a large uncertainty in the physics involved in the formation and survival of dark 
matter substructures. 

In Fig. ^we plot the flux of 7-rays in units of (GeV cm^ s sr)^^ arriving on 
Earth averaged in all directions for =100 GeV and with {av)Ann = 3 x 10^^^ 
cm'^/s. The curves refer to the 7-ray flux due to the full dark matter profile, made 
of the smooth and clumped components. For the type II scenario, the 7-ray flux 
contributed by the clumped component is comparable to the contribution of the 
smooth dark matter profile, while for the type I scenario, the clumpy component is 
overwhelmingly larger than that due to the smooth component. 

The fluxes plotted in Fig. ^ are obtained choosing the minimum clump mass of 
Mc^min = IO^Mq, but the dependence of these fluxes on the value of Mc^min is only 
logarithmic. In both flgures the dotted, dashed and solid lines correspond to SIS, 
Moore and NFW clump density proflles respectively. 

Also shown are the EGRET data (straight line) on the extragalactic diffuse 
7-ray background which can be fltted from 30 MeV to ^ 30 GeV by 1^ dUdE ~ 
1.36 X 10^® (cfv) ^^'^ GeV^^cm^^s^^sr^^. Depending on the density proflle the 
fluxes have different scalings with the neutralino parameters and {<7v)Ann'- 

^SIS OC {<7v)](^^m^^/^ ^NFW,Moore OC (av) Annm'^ , (11) 

these scalings are the same for the type I and II scenarios. 

It is clear from Fig. ^that the comparison between our predictions and the ob- 
served diffuse background is meaningful only for the type I scenario, with highly con- 
centrated clumps. In the second scenario, the fluxes are too low, with the exception 
of the case in which the density proflle is the SIS one. For the other cases the region 
of parameters that can be constrained is already ruled out from accelerator exper- 
iments > 50 GeV -^Si and from theoretical arguments {crv)Ann < 10~^^cm'^/s 
1^ 
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Figure 2. Limits on the neutralino parameter space in the first scenario for clump density profile. 

The situation is different for the type I scenario. As anticipated above, aU the 
density profiles imply 7-ray fluxes comparable to or largely in excess of the EGRET 
observations of the diffuse 7-ray background. These observations may therefore 
be used as a tool to extract severe constraints on the neutralino parameter space, 
which unfortunately are restricted to the type I scenario. Using the scalings given 
above (cfr. Eq. we can determine the regions of the parameter space which 

are ruled out by our calculations. 

Using EGRET data at 10 GeV, we find that SIS clumps in our halo are ruled out 
in the region {(Jv)Ann > 2 x l{)-'^^{m^/imGe\ffcnr'/s for hQGeV < < 3TeV, 
and {av)Ann > 6 x 10-4°(mx/100GeV)^cm3/s cm^/s for > 3 TeV. 

Moore et al. clumps are also strongly constrained: the 10 GeV EGRET data 
require that {av)Ann < 2.5 x 10-28(TOx/100GeV)3/2cm3/s for bOGeV < < 
STeF, and {av)Ann < 10-30(TO^/i00GeV)3cmVs for m^>3 TeV. 

If we extrapolate the EGRET measurement of the extragalactic diffuse 7- 
ray background to 100 GeV, the bounds get tighter: for between 50 GeV 
and 30 TeV, the flux from clumps is below the EGRET data if {<Jv)Ann < 
6 X 10-29(TOx/100GeV)3/2 cm^/s, while for > 30 TeV, the region {av)Ann < 
10-32(mx/100GeV)3 cm^/s is allowed. 

The NEW clumps in the type I scenario are the ones that are more weakly 
constrained. The bounds that can be placed by EGRET at 10 GeV are as follows: 
if is between 50 GeV and 3 TeV, the allowed region is defined by {av)Ann < 
10-26(mx/100GeV)3/2 cm^/s. For > 3 TeV, the allowed region is instead 
{(Tv)Ann < 6 X 10~29 (m^/ lOOGeV)^ cm^/s. If we extrapolate EGRET data up to 
100 GeV, the bounds become as follows: for between 50 GeV and 30 TeV, one 
must have {av)Ann < 3 x 10''^'^ (mx/lOOGeV)^^^ cm^/s while for > 30 TeV, 
the allowed region becomes {<Jv)Ann < 5 x 10~^^(m^/100GeV)'^ cm^^/s. All these 
bounds are shown in Fig. |21 The scatter plot reported in Fig. |21 represents all the 
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admitted values of {m^, {av)Ann) allowed by SUSY models 



5. Dark Matter and Magnetic Field at the GC 

The possibility of detecting neutralino annihilation trough the 7-ray emission, as 
discussed in the previous sections, is not the only one. Another appealing possi- 
bility is related to the synchrotron emission of electron-positron pairs produced by 
neutralino annihilation at the GC. In this section we will discuss the major char- 
acteristics of the DM density profile at the GC, taking also into account the effects 
on the GC magnetic field of the advection flow produced by the central BH. 

The spatial distribution of dark matter in galactic halos is still a matter of much 
debate (see, e.g., 1^^. Numerical simulations suggest that coUisionless dark matter 

forms cuspy halos while some observations argue for a flat inner density profile!^. 

on 

The standard numerical dark matter halo is the NFW profile " which is expected 
to be universal. However, recent simulations have found more cuspy halos^^as well 
as shallower profilesl^SI At present, it is not clear if there is a universal dark matter 
halo profile, but the NFW profile seems to represent well the range of possibilities. 
Therefore, we assume that the NFW profile describes well the dark matter in our 
Galaxy and we will use the NFW profile with the parameters introduced in §2. 

There is now growing evidence for the presence of a supermassive BH at the 
GC, with mass ~ 2 x 10^ Mq. In fact most galaxies seem to have central black 
holes with comparable or even larger masses. The presence of a BH can steepen 
the density profile of dark matter by transforming the cusp at the GC into a spike 
of dark matter El. xhe density profile of the spike region, where the gravitational 
potential is dominated by the BH is described by 

Here, S is the slope of the density profile of dark matter in the inner region {S — 1 
for a NFW profile), and 6sp = (9 — 2(5)/ (4— S). The coefficients as and g{r) can be 
calculated numerically as explained in detail inB. It is possible to identify a spike 
radius Rsp where the spike density profile given in Eq. (|12|l matches the NFW 
dark matter profile. In other words, at Rsp the gravitational potential is no longer 
dominated by the central BH. 

Neutralino annihilations affect the density profile in the spike by generating a 
flattening where the annihilation time becomes smaller than the age of the BH. This 
effect produces a constant neutralino density given by Pcore = (av)T'' tsH ' ^'^^^^ 
Ibh is the BH age. The effect of annihilations on the spike density profile can be 
written as 

/ N P'spi''') Pcore , , 

which accounts for the flattening in the central region. 
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Several dynamical effects may weaken or destroy the spike in the GC BEll 
depending on the history of formation of the central BH. If the spike is not formed 
or gets destroy, the central region should be described by the cuspy profile such as 
in the NFW case. Here we consider both cases and show that the observed emission 
is stronger than the predictions for a NFW cusp, while the spike generates signals 
well above the observations. 

In order to determine the synchrotron emission produced by the electrons that 
come from neutralino annihilations it is necessary to determine the magnetic field 
present in the GC region. In what follows we will determine the magnetic field 
strength assuming the equipartition between magnetic field energy and the kinetic 
pressure due to the accretion fiow inside the central BH. 

We model the accretion flow of gas onto the BH at the center of our Galaxy 
following a simple approach described inl^. More detailed models of the accretion 
flow around the BH lead to corrections which are negligible when compared to 
the uncertainties in the dark matter distribution. In the model we adopt, the BH 
accretes its fuel from a nearby molecular cloud, located at about 0.01 pc from the 
BH. The accretion is assumed to be spherically symmetric Bondi accretion with a 
rate of mass accretion of Al = 10^^M22 g s~^. The accretion onto the BH occurs 
with a velocity around the free-fall velocity, such that 

v{r) = ^2GMBH/r = c{-^\ (14) 

where Rg = 2GMbh/c^ = 7.4 x 10"(AfBH/2.5 x 10^ Af©) cm is the gravitational 
radius of the BH and Mbh is the BH mass. Therefore, 

/ n r \ 1/2 / \ -1/2 



v{r) = 1.46 X 10« cm s^^ (15) 

^ ' V2-5 X 106^/0/ VO-Olpc/ 

Mass conservation then gives the following density profile: 

M M f r ^ "^/^ 



i\/r \ -1/2 / \ -3/2 



such that 

p(r) = 5.6 X 10~^'M22 ( ) ( — '■ ) g cm"^ . (17) 

' V2-5xlO6M0y VO-Olpcy ^ ^ ' 

Following we assume that the magnetic field in the accretion flow achieves 
its equipartition value with the kinetic pressure, namely pv'^/2 — i?(r)^/87r. With 
this assumption. 

It is believed that magnetic fields in the accretion flow will in general reach the 
equipartition values described in Eq. (jl8|l . However, smaller fields may be reached 
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if the equipartition is prevented somehow. In what foUows we will always assume 
the equipartition field in deriving the synchrotron signal from electrons produced 
by neutralino annihilation. 

6. Neutralino Annihilation at the GC 

Following our discussion of section §3 we will assume that neutrino annihilation 
channel is dominated by quark-antiquark production, wi th a large production of 
pions described, as already discussed in §3, by the Hill 1^ spectrum. 

The spectrum of electrons (and positrons) from the tt* decays is calculated by 
convoluting the spectrum of pions and muons. For relativistic electrons the electron 
spectrum reads 

ml 1 dne{Ee,E^,E^) 



'^Jb™ " ' "'ml- ml - ml dE^ 

(19) 

where, neglecting the muon polarization, we get 



dneiEe,Ef„E„) _ 1 



dEe E^I3 



2 f 3+13^ 
3^ I 1-/3 



if 0<.<i^, 



(20) 



with e = Y^^- Here /3 is the pion speed and E^ and iJ^ arc the total energies 
of electrons and muons respectively. The two limits of integration E'^^'^[E^) and 
E^^^{Ef^) can be derived by inverting the following equations: 

E,<^[ml{l + f3) + ml{l-P)] (21) 

E,>^[mlil-f3)+ml{l + f3)] . (22) 

Finally, the injection of new electrons, produced in neutralino annihilation, 
at the distance r from the BH and at energy E can be written as Q{E,r) — 
{l/2){pi:)M{r)/m^)'^We{E){av)Ann, where the density of dark matter has the profile 
PDMir) discussed in §5 (NFW or spike). 

The spectrum of particles at a position r in the accretion flow around the BH 
is the result of the injection on newly produced electrons at the same position, 
radiative losses of these electrons and the adiabatic compression that may enhance 
their momentum while they move inward. Here we neglect spatial diffusion, which 
occurs on larger time scales. 

The transport equation including all these effects can be written as follows: 

- [rMr)] P^ + plQ-p [P'pir,p)f] = Qir,p), (23) 

where f{r,p) is the equilibrium distribution function of electrons injected according 
with Q{r^p), and losing energy radiatively as described by the function p{r,p) = 
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dp{r,p) / dt. Here v{r) = —cirjTlg)^^!'^ is the inflow velocity. The equation can be 
solved analytically if the electrons remain relativistic everywhere in th e flu id. The 
assumption of relativistic electrons can be safely used as discussed in In Eq. 
the term 

1 19 

Vad = -gPVw(r) = --^P^ [f^ w(r)] (24) 

describes the rate of change of momentum of a particle at the position r due to 
adiabatic compression in the accretion flow. The rate of adiabatic momentum en- 
hancement should be compared with the rate of losses due to synchrotron emission: 

Psyn[r,P) = -(JT 7 , (25) 

where ctt is the Thomson cross section and 7 is the Lorentz factor of the electron. 
The magnetic field Beq{r) depends on r as described in Eq. H18|) . 

In order to solve the transport equation, Eq. H23|) . we discuss different loss 
processes using the equipartition field Beq. We first consider synchrotron losses 
(jPsyn), followed by inverse Compton losses (pics) and synchrotron self-compton 
scattering (pssc)- The rate of synchrotron losses from Eq. (|25() is given by: 

Psyn{r,p) c = 1.6 X 10~^^ ( oipc ) 7^ ^""S ^^'^ Thomson regime, losses 

due to Inverse Compton Scattering (ICS) off a photon background with energy den- 
sity Uph has the following form pics = ^crxUphj'^ , ICS dominates synchrotron 
losses only if > SirUph- Assuming that Uph is independent of r (i.e. a fixed 
photon background) ICS becomes important at large radii (^ O.Olpc) and only if 
Uph > lO^eV cm~^. Such a strong photon background is unlikely to be present 
at the GC region. For comparison, the CMB radiation has Ucmb ~ 0.25eV cm~^ 
while the optical background has Uopt ~ leV cm^^. If ICS off a fixed background 
is not dominant at large radii, it becomes even less important as small radii when 
compared to synchrotron losses. Consequently, we safely neglect the role of ICS off 
photons of fixed photon backgrounds. 

The electrons, radiating in the strong magnetic field near the BH generate a 
photon background that can become quite intense. The rate of losses due to ICS 
of electrons off the photons generated through synchrotron emission by the same 
electrons is pssc(,f,p) = ^'^rU^fJ^ {r)j-^ , where the photon energy density generally 
depends on the radius r. The photon density Up^\r) is a nonlinear function of the 
distribution f(r,p). In other words, the term p in Eq. (|23|) depends in turn on f{r,p) 
when synchrotron self-Compton scattering is included. If this contribution cannot 
be neglected, an analytical solution of the transport equation becomes unattainable. 

Given the distribution function f(r,p), one can calculate the synchrotron emis- 
sivity jivjr) (energy per unit volume, per unit frequency, per unit time). The 
photon energy density at the position r is then proportional to the integration over 
all lines of sight of the emissivity, with the possible synchrotron self-absorption 
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Figure 3. Electron density per unit energy as a function of energy (solid lines) at r = 10^ Rg 
(upper curve) and r = 10^i?g (lower curve). Superimposed (dashed curves), we plot the function 
n{E, r) obtained without the effect of advection. The left panel illustrates the spike case, while 
the right panel applies to the NFW case. 



taken into account at each frequency. However, the distribution function f{r,p) is 
not known a priori, and the problem becomes intrinsically nonlinear. The approach 
that we follow here is start with neglecting synchrotron self-Compton scattering and 
check a posteriori whether the assumption is correct in the situation at hand. 

An analytic solution of the transport equation can be derived when p{r,p) is 
dominated by synchrotron losses as in Eq. (|25|l . In this case, the equation admits 
the following analytical solution: 

^'•'■'> ^ J {iYT" (it)' (If ■ w 

The function pinj — Pinj [Pi Rinj] corresponds to the injection momentum of an 
electron injected at the position Rinj that arrives at the position r with momentum 
p. This injection momentum can be obtained by inverting, with respect to pinj 
the solution of the equation of motion of the electron, in the presence of adiabatic 
compression and radiative losses: 

^^h(^y\2_^ (z_y\ (27) 

dr c \RgJ ^ 2Rg ^ \rJ ^ ' 

The solution of this equation, with initial condition p[r — Rinj, Pinj, Rinj] = Pinj is 



P\t, Pinj : Rinj] — Pinj 



2fco R^ 
3c 



~Pinj 



R 



in] 



3/2' 



R. 



in] 



1/2- 



(28) 



In the absence of synchrotron energy losses, particle momenta only change due to 
adiabatic compression, and the momentum of a particle changes according with the 
well known p = Pinj{r / Rinj)~^^^ , valid for the case of free fall. 

The joint effect of the energy gain due to the adiabatic compression and the 
energy losses due to synchrotron emission generates a new energy scale Pm in the 
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system: pm = 2koR^ ~ 3Trr{mc) ^ ^j^gj-g _ ar Bq / Gtt (mc)^ (ko — no synchrotron 
losses). Introducing p„i we can rewrite p[r,pi„j, i?i„j] as: 




3/2\ / ^ X 1/2 



Rinj J j \R 



(29) 



From this expression it is clear that, at any fixed position r, adiabatic com- 
pression dominates over synchrotron losses if the injection momentum is lower than 
Pm- In this case, the electron energy increases while the electron moves inward, 
until the rate of synchrotron losses become important. The opposite happens when 
the electrons are injected at momenta larger than p™, since synchrotron losses are 
important from the time of injection. The momentum pm can be interpreted as the 
momentum where the two competitive effects of adiabatic heating and synchrotron 
losses balance each other. Thus, particles accumulate at momentum pm- This phe- 
nomenon depends on the distance from the galactic center: at large distances from 
the BH the momentum pm, which scales linearly with radius, is large and the lo- 
cal rate of injected electrons is low, therefore, the accumulation is small. At small 
distances the accumulations at Pm grows. 

We can define the electron equilibrium spectrum n{E, r) which is related to 
/(i?, r) through the relation n{E, r)dE — ATrp^f{p, r)dp. In Fig. |31we plot n{E, r), 
in the two cases of an NFW (right panel) and a spike (left panel) density profiles, 
as a function of the electron energy at two different radii, r — lO'^Rg (upper solid 
curve) and r — 10*i?g (lower solid curve); we have fixed = 100 GeV and 
{crv)Ann = 10^^''cni^/s. The dashed lines illustrate the solution of the transport 
equation when adiabatic compression is switched off and only synchrotron losses 
are included. The accumulation effect described above manifests itself through the 
appearance of the spiky structure at momentum pm- The accumulation is less 
pronounced at large radii, as expected. 

We conclude this section by addressing the issue of the synchrotron self- Comp- 
ton scattering. As explained above, this effect cannot be accounted for in an analyt- 
ical approach to the transport equation, since it is intrinsically nonlinear. Instead, 
we check a posteriori if neglecting SSC was a good assumption. The photon energy 

arguments: 



density as a function of r, U^fJ^r), is easily calculated on the basis of symmetry 



u;rir) = - 



dr'r" I di'j{u,r )~{- / dr / di'j{u,r ) 



oo 



(30) 



where j{v,r) is the synchrotron emissivity. 

This energy density can now be compared with the magnetic energy density 
at the same location, _B^(r )/87r. The results are plotted in Fig. 01 the solid line 
represents the photon energy density Upf"{r), while the magnetic energy density is 
plotted as a dashed line. The calculations are carried out for a dark matter density 
profile with the spike at the center and a neutralino mass of 100 GeV. The magnetic 
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Figure 4. Comparison between the energy density of the synchrotron emitted photons f/p^ (sohd 
Une) and the energy density associated to the magnetic field B^/Svr (dashed line). 

energy always dominates over the photon energy, although at large distances from 
the BH the difference between the two curves reduces to about one order of magni- 
tude. The curves in Fig0]are obtained without taking into account the synchrotron 
self-absorption effect, therefore the calculated photon energy density should be con- 
sidered as an upper limit. Thus, neglecting SSC is a good approximation for the 
present scenario. 

Let us conclude this section addressing the issue of the synchrotron self- 
absorption (SSA) mechanism. Synchrotron radiation can be reabsorbed by the 
radiating electrons when the system is sufficiently compact. This phenomenon is 
particularly effective at low frequencie s. 

Following the standard procedure 1^ we have included the SSA effect in the 
emission evaluation finding that this effect is efficient onl y in the case of the spike 
density profile and only for frequencies v < 10^^ Hz This results holds only 
in the case in which both synchrotron losses and adiabatic compression are taken 
into account. On the other hand if only synchrotron losses are considered the 
SSA effect becomes important even at lower frequencies. The dependence of the 
SSA on the neutralino parameters scales with n{E,r) cx ((Tu)yi„„mx'^^^. At fixed 
{<Jv)Ann, the electron density drops with increased m^, and the SSA effect decreases 
as well. Finally, in the case of an NFW density profile without a spike, the rate of 
electrons injection remains such that the effect of SSA can always be neglected at 
the frequencies of interest. 

7. Synchrotron Emission from the GC 

In this section, we present our results in terms of the synchrotron luminosity pro- 
duced by relativistic electrons from neutralino annihilations in the GC. We have 
considered both cases of a spike density profile and an NFW density profile. Taking 
into account the effect of SSA, that is efficient only in the spike case as discussed 
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before, we have computed, through the equihbrium spectrum of electrons n{E,r), 
the synchrotron luminosity integrated over the all sky (i.e. over all lines of sight). 

In Fig. El we fix = 10^ GeV and {(Jv)Ann = 10"^"^ cm'^/s and show the 
luminosity obtained with both the spike (left panel) and NFW (right panel) density 
profile for two cases: the continuous line shows the case in which advcction is taken 
into account, while the dashed line shows the case in which only the synchrotron 
energy losses are present. Comparison of the two curves in Fig. El shows that the 
resonant behavior due to the combined effect of advection and synchrotron losses 
becomes important for frequencies up to 10^"* Hz. In this frequency range, the 
electron density is much higher when advection is included as compared to the 
pure synchrotron case. This effect produces an increase in the emitted luminosity 
of about one order of magnitude. In addition, we can see that SSA decreases the 
emitted radiation in the frequency range 10^° — 10^"'^ Hz. Moreover, the SSA effect 
becomes relevant only when advection is taken into account and only in the case of 
the spike density profile. 

In order to constrain the shape of the dark matter density profile at the GC or the 
parameters of the dark matter particle, we compare the calculated luminosities with 
the observations. The well-studied source at the GC, Sgr A*, has been observed in a 
large range of frequencie s, fr om the radio up to the near-infrared. Experimental data 
from 10^ Hz to 10^^ Hz I""' are displayed in Fig. |B1 In the same figure we show the 
two cases of spiky (left panel) and NFW (right panel) density profile, for = 10^ 
GeV (upper curves) and — IQ^ GeV (lower curves) and (crw)^„„ = 10^^^ 
cm'^/s. In Fig. Owe also show the results for the cases in which only synchrotron 
losses are taken into account (dashed curves). From Fig. (HI it is clear that a 
spike in the GC induces much stronger emission than the observed flux. Therefore, 
either there is no spike in the dark matter profile or neutralinos are not the dark 
matter. This conclusion agrees with previous studies ^ that show how changes 
to the neutralino parameter do not gap this large discrepancy. The parameter 
space for neutralinos m^, (crw)^„„ has been studied extensively!^. From accelerator 
studies = 10^ GeV is a lower limit and from cosmological constraints (i.e., dark 
matter density) and {av)Ann must lie in the range: € [lO^GeV, lO'^GeV] 
and {av)Ann G [IQ-^^cm^/s, IQ-^ecniVs] ■ 

The emitted luminosity scales with the neutralino mass and the annihilation 
cross section through the electron density n{E^ r), following the same relations ob- 
tained in Eq. Hll|l . This implies that the luminosity can only be lowered by at most 
one order of magnitude, by increasing the neutralino mass up to ~ I TeV. While it 
is easy to increase the luminosity by an order of magnitude without violating any 
bound simply by increasing the cross section of neutralino annihilation, it is hard 
to decrease it. In fact, as {<Jv)Ann is lowered well below 10^^^ cm^^/s, the density 
of neutralinos becomes smaller than the necessary density to account for cold dark 
matter. 

The large enhancement of the neutralino density due to the spike profile produces 
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Figure 5. Emitted luminosty in two cases: 1) advection and synchrotron losses (solid line); 2) 
only synchrotron losses (dashed line). The computation is performed with m^^ = 10^ GeV and 
{uv)^„„ = 10~^^ cm^/s. The left panel represents the Spike case while right panel the NFW case. 

a synchrotron luminosity that is difficult to reconcile with observations. Moreover, 
this conclusion only gets stronger when advection is included. The SSA effect does 
not affect the discrepancy between a neutralino spike and the observations, since 
it changes the luminosity only at frequencies in the range 10^^ 10^^ Hz. The 
situation is completely different for the case of the NFW density profile. In this 
case, the synchrotron luminosity is always less than the experimental data as can 
be seen in Fig. El NFW is consistent with observations even if {av)An7i = 10^^® 
cm^/s is considered. 

In this paper, we only considered two possibilities for the dark matter density 
profile: the less concentrated hypothesis (NFW) and the most concentrated one 
(spike). Other proposed profiles such as the Moore et al. profile should generate a 
luminosity in between the spike and the NFW cases. In this case, observations are 
likely to place more stringent limits on the neutralino parameters instead on the 
density profile itself. 

8. Conclusions 

In this paper we have discussed two possible avenues that lead to an indirect de- 
tection of the annihilation of neutralino CDM. We have discussed the effect of this 
annihilation in the halo of our galaxy as well as in its center. 

The expected emission from the halo annihilation drops into the 7-ray frequency 
range and it generates a diffuse 7-ray background. The Dark matter sub-structures 
of the galactic halo can be a dominant component of the diffuse 7-ray background 
depending on the concentration and location of DM clumps in the inner regions of 
the galactic halo. In order to bracket the range of possible fiuxes we have considered 
two extreme scenarios, that we named type I and type II. The first corresponds to 
extremely concentrated clumps present everywhere in the Galaxy halo, while in the 
second scenario the clumps are much less concentrated and are completely destroyed 
by tidal effects in the inner 10 kpc of the Galaxy. While the type I scenario allows 
one to put very strong constraints on the properties of neutralinos and on the density 
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Figure 6. Luminosity, compared with experimental data from Sgr A*, in tfie case of the spiky 
density profile (left panel) and NFW density profile (right panel). Dashed curves represent the 
luminosity obtained neglecting the effect of advection. The computation is performed with rrtx = 
lO'^ GeV (upper curves) and mx = 10^ GeV (lower curves). The annihilation cross section is 

(o-f>Ann = 10-2''' cm3/s. 

profile inside clumps, the type II scenario generates fluxes of diffuse 7-rays which 
are barely detectable, with the exception of the case in which the density profile 
of the clumps is modeled as a SIS sphere. For the type I scenario, most of the 
parameter space of neutralino dark matter is ruled out if the density profile of dark 
matter clumps is in the form of a SIS sphere or a Moore profile. Weaker bounds can 
be imposed on the neutralino parameter space in the case of NFW density profile. 

The second possibility that we have considered here is related to the electron- 
positron emission in the GC magnetic field. Electron-positron pairs are expected 
to be copiously produced in neutralino annihilation, from the comparative study of 
their synchrotron emission at the GC with the experimental data, already available 
in the radio and near IR frequency range, we have obtained severe limits on the 
neutralino density profile at the GG. Using our results we can reaffirm that a spike 
density profile is ruled out by these radio and near IR observations. We reached this 
conclusion by a careful consideration of the accretion flow and the loss processes 
in the transport equations for the neutralino generated electrons and positrons as 
well as the radiative transfer. However, rejecting the spike hypothesis the observed 
emission from the GC is consistent with neutralinos following an NFW density 
profile. We can conclude stating that other proposed density profiles, like the Moore 
et al., should generate an emission in between the two profiles we have considered 
here. In this case the experimental data from Sg A* will put severe constraints on 
the neutralino parameters. 
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